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ABSTRACT 




The coarsening behaviour of carbonitnde precipitate 
has been examined in a high strength low alloy steel (HSLA) 
containing 0 105 wt / Ti and 0 089 wt/ V as carbonitride 
forming elements in the temperature range 733-770*^0 The 

steel was aged for periods ranging from 1 h to 343 h at 
733°C» 752°C and 770^^0 temperatures Transmission Electron 
Microscopy was carried out on carbon extraction replica 
Particle size distributions (PSDs) were determined for 
various aging times at the three aging temperatures 
Although a^= kt rate Law was observed during coarsening 
None of the presently available coarsening theories showed 
complete agreement with the experimental particle size 
distribution A normal distribution rather than a 

log-normal distribution was observed The value of volume 
diffusivity, Dv IS observed to be comparable with the 
corresponding published value 
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CHAPTER 1 

INTRODUCTION 

The mechanical properties of HSLA steel are strongly 
dependent upon the size and distribution of the carbide 

precipitate Since the carbide particle can coarsen during the 
initial heat treatment and during subsequent service^ tt is 
important to be able to predict the Kinetics of coarsening The 

ideas of precipitate coarsening was first formulated by Qstwald 
C13 as early as in 1900 The precipitate coarsening describes the 
change in average particle size with increasing time It results 
from higher solubility of small particles, which may tend to 
dissolve in the matrix and cause the growth of larger particles 
The driving force is provided by the reduction in total 

xnterfacial energy ot the system when the average particle size 
increases with time This can be best illustrated with the help 
of the free energy-composition diagram, Figure 1 1 Let us 

consider a large spherical particle of radius r being surrounded 
by a number of smaller particles whose agerage radius is r, Figure 
12 A common tangent construction between and Gp phases of 
the free energy - composition diagram established equilibrium 
between the two phases Such construction indicates that large 

particles have lower solubility of the solute in the matrix in 
which they are contained It is clear from the figure I 1 that 
when large particles grow at the expense of smaller particles, 
there is a reduction in the interfacial energy per unit volume of 
the material The energy thus released acts as driving force for 
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the coarsening reaction 

Because of the importance of the precipitate stability in 
commercial alloys the coarsening of precipitates has been studied 
experimentally by numerous investigators Detailed theoretical 
analysis of the kinetics of coarsening was carried out by 
Lifschitz and Slyozov C21 and independently by Wagner C33 in what 
has become to be known as the LSW theory of coarsening The 
principal finding of the LSW theory is that the average particle 
size increases with the cube root of the time of coarsening 

The solubility of a particle in the matrix as a function of 
its radius is best illustrated by the Gibbs - Thomson equation, 
written as 


in 

where 


m * 
RT 


dV 


(1 I > 


r a/ (i 
CD ex/ ft 


m 

r 

dA 

dV 


For spher 
A = 4rTr^, 
So the term 
substitution y 


Solubility of a particle of radius r, 

Solubility of a particle of infinite radius, 

Molar Volume of the precipitate. 

Interfacial free energy, 

Rate of change of surface area of the precipitate 
with volume and R & T have their usual meaning 
ical particles of radius of h 
V = 4/3 irir^ 

dA/dV in equation (1 1) equals 2/r which upon 
ields 


r^a/ft 


(D^CA/ft 


') 


2V^^ 

m' 


In 


RTr 


(1 2 ) 


1 


Assuming typical values of the terms an the right hand side 
of equation <1 2) 

Y = 100 mJ/m^ 

R B 314J/male'^K 

T - iOOO'^K 

to X 10 ^m^/mole 

m 

r = 10 ^ ta 10 


Then 




2 X too X 10 X 10 ^ K 10 ^ = 2 4 X lo"^ 
8 314 X 1000 X 10“^ 


Since the term on the right hand side has a small value, 
equation <1 2) is written as 


1 r a/ 0 _ oSyOt//?! _ 

J " 


2 A 

m' B 


<1 3> 


If an isolated particle of the precipitate phase is rich in 


solute then the concentration distance profile in the early stages 


of its growth is as shown in Figure 1 3 However, nucleation of 
precipitate particles in the vicinity of an isolated particles and 
their continued growth make the diffusion field overlap The 
particles stop growing thereafter New region cannot nucleate 
solute-rich ft phase because the neighbouring matrix is depleeted 
of solute The smaller particles dissolve during coarsening and 
the solute thus released moves down the concentration gradient 
(built gradually from smaller to larger particles) to attach 
themselves to larger particles The proces^of coarsening 
therefore is one where some particles are continuously dissolving, 
some are continuously growing and some remain quasi-stat lonary at 
any given small time interval During this interval some 



Matr X 



1 3 Concentration - distance 
which IS growing 


profile on 


a particle 



1 4 Theoretical particle size distribution as 
predicted by the LSWC5], BWC8], BWEMC93, and 
V6C10] theories The steady-state distribution 
are presented for the BW and BWEM theories at 
precipitate volume fraction of 60/ and for the V6 
theory at 50 Vol / CRef 153 
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particles must have dissolved completely and disappeared The 

particles which remain quasi-stationary , i e ^which is neither 
growing nor dissolving at any given small time interval is termed 
as the critical size particle of radius r It is understood that 
particles that remained quasi-stat lonary in a given interval may 
not be so in subsequent intervals This is so because with the 
coarsening of particles^ the average particle size increases and 
the critical particle size also increases with it The net 

results being that the number of particles remaining after any 
time interval decreases with time and the average particle size 
increases with time 

There are a number of coarsening mechanism These are 

1 Volume diffusion controlled 

2 Ledge mechanism 

3 Grain boundary diffusion - controlled 

4 Diffusion through dislocations 

5 Interface controlled 

When the rate controlling mechanism is the diffusion of 
solute through the regular lattice sites^ the coarsening mechanism 
IS said to be volume diffusion controlled (VDC) In order to 
establish a functional relationship between the particle size and 
the time of coarsening^ It is essential to make some simplified 
assumption These are 


1 

The concent rat x on 

grad i en t 

IS considered 

to 

be 

time 


independen t 

This 

will allow 

us to use the 

Pick s 

f 1 rs t 

Law 


to adequately 

desc 

ribes the 

concentration 

grad i 

ent 

under 


steady state conditions 
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the a./ ft interface 

^aft 

of curvature 

This 

1 S ( 

of lOoS or so 

Since 

the 


IS 


2 The salute diffusion coefficient* D is independent of 
concentration of the solute This condition will again be 
fulfilled because dilute solutions are generally used in 
precipitation studies 

3 The interfacial energy 

independent of the radiu? 

critical below particle si 

of the particle at the termination of growth is of the order 

of lOOS* the interfacial energy will not be seriously 

affected during coarsening Thus a constant value of i> ^ 

afi 

will be assumed for the entire coarsening process 

4 For spherical particles, the effective diffusion distance is 
considered to be proportional to the radius of curvature of 
the particle The value of the proportionality constant, <m, 
which iap> reality is a complicated non-analy tical function of 
several variables is assumed to remain constant during the 
entire coarsening process 

^ Diffusion Controlled Coarsening 

For spherical particles of radius r, the flux to the 

particle fflwom the matrix is given by the Fick s first Law as 


= - D 4fTr 


2 ^ 

Ax 


(1 4) 


Where D = rate of diffusion of solute through the lattice 
2 

4fTr = surface area of the sphere 
Ax = effective diffusion distance 

From the Gibb 'Si Thomson s effect, equation (I 3), the solute 
concedntrat ion difference AXg is given as 



and 


i 


AX 


B 



r ^ ^ 

r a/ fi r a/ ft 

^ m 

i ^ i 


B " 

' B RT 

r — 


^ j 


r 


Ax effective = ar 

which upon substitution in equation (1 4) yields 


2vV^ 


<n 


D 4nr^ “x"'''’ - , 

® OiRTr^ 




(1 5) 


The flux can be equated to the flux arising out of the 
growth of spherical shaped particle from the matrix which gives 



- 

dV 

^X^ - 

■ x“l 


B 

d t 

B 

"J 


. 2 d r _ 

at ■ ■’i 


<1 i* 


Where 


V = Volume of the sphere and 


[xb - 


* Solute concentration difference across the ct/ft 
boundary 


, Upon equating the flux equation for and and rearranging 
terms^ we gat 


dr 

dt 


2DV^ i> 

m B 


o<RT <X 


B 


X^)r^ 


1 - H 


<1 7) 


On rearranging the terms of equation (1 7) and integrating on 
both sides we get 




2D V "Xg 


V m 


<l-£) 

a 


0 "'’T <*B - 


dt 


<l 8) 


Treating as constant and applying results of the LSW theoryj 
equation (1 8) becomes 


_ _ _ ®x2^" 

—3—3 8 m B 

r - r = « D 

o 9 V ^ct 


(1 9) 


ctRT ( X 


B 
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Where r is the average particle radius after time t, r is 
the average particle radius at the onset of coarsening is the 
diffusion coefficient of the solute in the matrix and other terms 
have already been defined 

D Ramkrishna and S P Gupta C43 derived the equation for 
platelets shape particle by volume diffusion controlled 
coarsenning They consider square platelets of edge length a and 
thickness h If the interfacial energies of the broad face and 
the edge with the a matrix are and respectively, application 
of the Gibbs Theomson equation yields 


.n = 

Ko J 


V' 4hv + 4aVc- 
m e F 

RT 2 ah 


(1 10 ) 


Where Xg and are the solubilities of particle of any size 
and a very large size respectively and is the molar volume of 
the ft phase 

From the Gibbs-Wulff theorem, The following can be written 


h/2 a/2 


( 1 11 ) 


Where W is the Wulff constant, which on substitution in 
equation <1 11) gives the solubility difference between particles 


of edge length a and mean value a as 
4 V X f, , I 

° V — [^ ; ) 


(1 12 ) 


The flux to the particle from the Fick s first Law is 


2 ^e'^m 

^ <2a ■^4ah) ^ta(a 7 2T a-- 


(■-tl 


<1 13) 


Where aa/2 is the diffusion distance, is the volume 

2 

diffusion coefficient and 2a + 4ah is the area through which 


The flux IS equated to the flux 


diffusion takes place 



arising out of the growth of particles as follows 


J 


2 





, “X“ 

■D <2a^*4ah) — 1-^ ^ 

V 2 

RTcta /2 


1 - 


(1 14) 


Letting h = a//9 > differentiating and rearranging yields 
a/2 ^ t 


1 


(a/2> d (a/2) 


a /2 


( 1 - ^ ) 


..ft ca^a 

= ~ I D +2) dt 

J 3 V ‘ RTot 


(1 15) 


o a 

The parametei ft is related to the specific interfacial 
energies of the two interfaces and can be calculated Treating oi 
as a constant and applying results of the Lifshi tz-Slyozov-Wagner 
theory gives 




■+2) D 

V 


e m B 


otRT 




(X'3-X«) 


<1 16) 


Except for the numerical factor 8/27 (.ft +2), equation <1 16) 
is identical in form with the kinetics equation of Wagner <3) 
Substitution for ft - l(h=a> enables the kinetics equation for the 
coarsening of cube-shaped precipitates to be obtained 


a 

2 


' H ■ 


a V g m B 
aRT <X^-Xg) 


(1 17) 


(a) -(a^) = ^ 


,, D V 
64 V e m B 


«RT <X^-Xg) 


<1 18) 


Equation <1 18) can be written in an approximate form 
a = (Kt)^'"^ 


(I 19) 



n 


where 


K<Reite constant) = 


D V 

64 V g m B 
^ ctRT (X^-X") 


1 2 Partlcal Size Distribution 

Given sufficient time of coarsening of the precipitate after 
initial growth, particle which are larger than the average will 
grow and those smaller than the average size will gradually shrink 
and disappear Thus a particle size distribution function, f<r,t) 
can be defined that corresponds to the number of particles of a 


given radius r at time t The total number of particle in the 
system at any instant is given by 


"p ^ 1 


f < r, t )dr 


initially, n^ is constant and equals the number of particles at 


the end of the growth process of the precipitate 


However, as 


some particles grow, some dissolve and others remain stationary 
during coarsening, the particle size distribution broadens The 
distribution of particle sizes evolves towards a quasi-steady 
state distribution, f<r,t), which is independent of the initial 
distribution f(r,0) The distribution function is written as 


f<r,t) = f(t) p h (p) 


(1 20 ) 


Where f (t) is function of time only 


and h(p) as solved by Wagner for a Gaussian distribution is 


h (p) 


W 


f 4 

exp 


3/2 - p 


<1 21 ) 



iz 


From which 

h{p) =0 for = 3/2 

p bexna one cut-off radius ratio also known as the maKimum 

2 

relative radius of the particle Taking f(r,t)a p h<p>, it is 
clear that, 

for p=0, f(r,t)=0 

and for p = 3/2 f<r,t) = 0 

This implies a sharp cut-off in th distribution, such that no 
precipitate particles exists with radius more than 1 5 times the 
mean radius of the particles The distribution has a maximum 
value at p = 1 135 

Fallowing Wagner s approach C31 , the integration can be 
carried out for 

P = ^ 2 yi®l-cl»'ng 

r r 



2 . 

r dr 



2D 


afi m B 


otRT 


( X^- 


x», 


or 


<r) - 


< r > ' 
o 


8 

9 


D -o 
V atft 


V^“’Xn 
tt\ 0 


«RT 


< 


x“l 


(1 22 ) 


<1 23) 


Although the cube rate law has been exhibited by a number of 
alloys The experimental particle size distribution IPSD s) do 
not concide with the prediction of the LSW C53 theory This 
discrepancy is mainly attributed to the fact that the LSW theory 
IS strickly applicable when the precipitate volume fraction is 
small and approaches zero obviously this zero volume fraction 



o 


approximation is not valid for many material systems As a 

result, considerable effort has gone into modifying the basic 

theory For example, L i f sh i z-Slyozou Encounter Modified <LSEM) 
theory It has been shown by Davies et al C71 which considers 
encounter between precipitate particles, the Brai Isford-Wynblatt 
(BW) E83, the Brai 1 isf ord-Wynblatt Encounter Modified (BWEM) C71, 
and the Voorhees— G1 icksman (V6) CIO] theories 

The theoretical PSD s predicted by the LSW C5], the BW CBl 

BWEM C9], and the VG CIO], theories are shown in Figure I 4 The 

latter three theories have attempted to incorporate the influence 
of finite precipitate volume fraction Figure 1 4 shows the 

steady state distribution which are predicted by the BW and BWEM 
theories for alloys with 60 7 qI / precipitate and by the V6 theory 
for allays with 50 Vol/ precipitate In the present study, the 
experimental PSD s were determined as a function of time and 
temperature 

1 3 Effect of Volume Fraction of the Prclpltale 

The theory of particle coarsening by Lifshitz and Slgozov and 
Wagner (LSW), is applicable only when the volume fraction of the 
precipitate is very small, under the condition, the mean distance 
between particles is large compared to the dimension of the 
particle Since the kinetics of coarsening is controlled by the 
diffusion of salute from a particle which is dissolving to a 
particle which is growing, it is expected that the mean distance 
between particles which is controlled in turn by the volume 
fraction, of the precipitate will influence the rate of growth 
As increases, the mean separation between particle* decreases 
and the distance that solute atoms have to diffuse on an average 
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to become a part of the growing precipitate become shorter A 

detailed analysis of the problem has been treated by Ardell C51 by 
considering a distribution of particle sizes surrounding an ith 
particle of radius r^, as shown in Figure 15 If the mean free 
path between the ith particle and its neighbour is X then, for a 

poly dispersed assembly of particles, the distance r at which the 

salute concentration in the matrix approaches the salute 
concentration of the particle - matrix interface of average size 
particle is given by 

r = r -*■ J (1 24) 

In order to introduce the effect of the volume fraction on 

the rate of coarsening, a relation between r and must be 

obtained This relation is difficult to obtain for a polydisparse 
assembly However, when the particles are considered to be all of 
the same radius r , the center to center distance, d, between a 
particle and its nearest neighbour is 

f -I 

' ’ ^ r * ^ '"n 


where the gamma function. 


CO 


r(V 


8v, 


-2/3 


e dy 


Now, from figure 15, X = d - 2r 

BV^ 


3^yv 


<V^>r 


amd r =» r + 


* f 

r e 


2 


a 26 ) 


T(V^) 


(1 27) 





I6 


r _ 1 ^ 1 
or - - 1 + ^ 


(1 2S) 


where n = 


2 3^yv^ 

e®'^f <v^) 


and P =' 


The effect of volume fraction on the parameter y) is shown in 
Figure 1 6 

Rate equation which is as fallows 


3 6D i, “x" 

dr _ m B 


dt 


RT <X' 


X") 


<l + 7?P) <p-l) 


or, 


V ^ ®X2 

_ a(i m B 

a p,— . w/? _ y^l 


(1 29) 


RT (X' 


B 




Where r is the critical radius of the particle in the initial 
o 

distribution 

r IS critical size which is neither growing nor dissolving 


V = 


P = 


• n ,,fi ®„a 

6D V . V‘ Xq .. 
of? m B dt 


RT <xj^ - Xg) dr*^ 

— which upon substitution in equation (1 14) 


,r, ijfi <»v^ 

_ _ 6D V - V' X-s T 

, — .3 /— v3 m B , — .3. 

<r) - (r ) * — ip) t 

° RT <X'q - Xg)^ 


<1 30) 


The value of has been calculated by Ardell C53 and it is 

found to have weak dependence on the volume fraction of the 
precipitate In the range that most studied of coarsening are 
carried out, p increases from approximately 0 94 to 10 as the 
volume fraction decrease from O I to zero 
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1 4 Estimation of O amd 

Equation (I 19) is of the form y = where indicate 

— 1/3 

the slope of line From Plot of a vs t atdifferent temperatures 


and Dv can be calculated The Arrhenius equation is 


D = D exp 

V Q ^ 


-q/RT 


a 31) 


where is the diffusion rate, is the rate constant 
Q IS the activation energy and 
RStT have their usual meaning 


Upon substitution of the value of D in f»quation (I 19) we get 


, ,/? oo^c*. 

^ ^m S 

otRT (X^ - Xg) 


-Q/ RT 


(1 32) 


This can beitivritten as 

, . D l> V' X_ o/o-r 

i-T 64 o m B -Q/RT 

KT = ^ ^ exp 

oR (X'q - X“) 


(1 33) 


or KT = A exp 


-Q/Rl 


or In KT = In A 


where 

, ,/9 oo^a. 

__ M °o‘' '4, ’‘b 

’ od (X^ - X“) 

We can determine the diffusion parameters D and Q from the 

Q 

intercept and ^lope of the straight line respectively from the 

Arrhenius plot On compalring the theoretical and experimental 

values of D and Q one can predict the mechanism that is operating 
o 

for the growth of precipitates 
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CHAPTER 2 

EXPERIMENTAL PROCEDURE 


2 1 Material 

The present in v'est igafc lan was carried out to stud/ the 

kinetics of precipitate coarsening in a HSLA Steel containing Ti 

and V The HSLA steel containing 0 1/ Vanadium and 0 08/ Titanium 

was obtained in the hot forged condition from Climax Molybdenum 

Corporation USA Table 2 1 shows the composition of this steel 

in the as ~ received condition It contains small amount of 

nitrogen which makes Vanadium nitride and Titanium nitrite 

precipitation possible during isothermal aging The elements Si 

and A1 are in very minor quantity and thus will have negligible 

influence The sole function of these two elements will be in 

contributing to solid solution strengthening The addition of Mn 

retards the rate of growth of ferrite from Austenite 

The steel had been hot forged to obtain rounds 18 mm in 

diameter and then cold rolled with a number of intermediate 

annealing treatments to obtain section 3 mm thick Small square 

2 

specimen size 6x6mm were solution treated for 9 minutes 
Isothermal transformation was carried out at temperatures 733'’c, 
752°C and 770°C in sealed quartz tube under vacuum for 1 h to 15 
days 

The experimental procedure and the techniques used for the 
study of kinetics of precipitate coarsening are as described in 
the following sections 



Table a 1 



Composition 

of HSLA 

Steel Cin 

Wt JO 


Element 

C 

N 

Si 

A1 

Mn V 

Ti 

Amoun t 

0 1 

0 016 

0 09 

0 048 

mmm 

in 

o 

o 



(a) 


(b) 




(c) 




Fiq 2 1 


Th» extraction replica technique 
surface <b) Etched to expose second 
Carbon coated, and <d) re-etched 
floated off CRef 113 


Ci) Polished 
phase, (c) 
and replica 




2 2 Replica Preparation 


In order to study the kinetics of precipitate coarsening of 
TiVCN the particle size and shapes were derived from carbon 
extraction replica technique This technique is illustrated 

schematically in Figure 2 1 The two phase structure to be 

replicated was first etched sufficiently in a nital solution 
containing 3/ nitric acid to expose the second phase A thin 

layer of carbon coating was done in the usual way Shadowing of 
carbon is to be avoided in order to ensure a strong continuous 
film and therefore, rotation of the specimen or diffuse deposition 
was done Subsequently shadowing was used to give information 
about the three dimensional shape of the extracted particles The 
specimen was next dipped in a nital solution containing 12/ nitric 
acid so that the matrix was attacked while the second phase, which 
was to be extracted remained unchanged The replica, containing 
particles of the second phase was collected on a copper grid The 
most critical part of this procedure was the second etching to 
free the replica and particles from the base metal It was 

necessary to score the carbon coating into squares before etching 
This floating carbon film containing particles was collected on 
capper grid (dia = 3mm) for further study 
2 3 Transmission Electron Microscopic Studies 

Carbon extraction replicas collected on copper grid were 
examined using the TEM < JEOL-2000FX ) In order to confirm the 
shape and size of the precipitate sufficient number of photographs 
were taken and in order to confirm precipitate phase is indeed 
Vanadium Titanium carbon i t r ide , selected area diffraction patterns 

were obtained from isolated large particle 
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2 4. Particle Size Measiir ement 

Edge length of the particle was measured manually with the 
help of a vernier caliper, having least count O 002cm From each 
samples measurements were made on approximately 500 particles 
2 S Actual Magnification Calibration 

The electron photomicrographs taken on 35 mm camera don t have 
magnification marker on them In order to carryout the 

calibration the same features were recorded on both the 35 mm 
camera and plate camera at a number of magnifications The plat© 

camera had the magnification marker on it The same precipitate 
particle size is measured in both the photograph and the actual 
magnification of 35 mm camera determined 
2 ® Particle Size Dtstribution 

The particle size distribution data were obtained by feeding 
the raw data to a computer programmed to give mean particle size, 
size distribution and particle size ratio, a/a for each specimen 
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CHAPTER 3 

RESULTS AND DISCUSSION 


3 1 Morphology 

A photomicrograph of the precipitate during coarsening is 
shown in Figure 3 1, which is a transmission electron micrograph 
of extraction replica derived from a specimen aged for 64 hours at 
770°C It clearly shows the square cross section of the 

precipitate In order to confirm it, the precipitate phase is 
indeed Vanadium-Titanium carbonitride (VTiCM), selected area 
diffraction pattern was obtained from isolated large particle as 
shown in Figure 32, as an illustration It indexes as Face 

centered cubic (FCC) lattice The cube geometry was maintained 
for the entire duration of coarsening as confirmed from the TEM 
micrograph of Figure 3 3 derived from a specimen aged for 216 
hours at 752*^0 During later stages of coarsening, there was a 
slight tendency for the edges of the cube to either get rounded or 
convert into small facets With these modifications, the cube 
shaped particles appeared as slightly rounded as shown in Figure 
3 4, derived from a specimen aged for 343 hours at 733^^0 The 

square geometry was not eliminated, however 
3 2 Coarsening Behaviour of Preclpitale* 

The coarsening behaviour of the VTiCN precipitate was 
examined as a function of temperature from 733°C to 770‘^C for the 
0 OQl wtY V, 0 0105 wt/ Ti, 0 IWt/ C and 0 0l6Wt/ N containing 
HSLA Steel The TEM photomicrographs of Figure 3 5 depict the 
change in the microstructure which occurs during aging at 770°C 



Fig 3 


Fig 



1 TEM Micrograph of extraction 
aged for 64 hours at 770®C 


rep 1 ica , 


specimen 



2 Electron diffraction pattern derived from 
isolated precipitate particle 





k- 0 bfim -H 


In TEM Micrograph Cube Shaped particles 
as slightly rounded 


appeared 





for times of up to 216 hours Figure 3 5<a) shows the initial 
micro structure after aging for 8 minutes The precipitate 
particles in this condition were cubical in shape and had an 
average particle size» a = 78 22 angstrom After aging for Ih, 
8h, 27hT 64h and 216h the precipitate particles increase in size, 
but some were coalescing as well Figure 3 5Cb) derived from a 
specimen aged for 216 hours at 770'^C shows the increased in 
average particle size 

Similar results were obtained during coarsening at 752°C and 
733'°C The precipitate particle wer^^ cubical in shape* as shown in 
Figures 3 6(a) and 3 6(b) derived from specimens aged for 8 hours 
and 64 hours respectively at 752*^C Photomicrographs presented in 
Figures 3 7(a) and 3 7(b) show the precipitate particle derived 
from a specimen aged for 8 minutes and 64 hours respectively at 
733°C 

3 3 Particle Size Ptstrlbutlons i 

The initial particle size distribution after 8 minutes aging 
at 770^C IS shown in Figure 3 3(a) The aging tim** was sufficient 
for complete transformation of austenite to ferrite at this 
temperature The austenite to ferrite transformation initiates 
after an incubation period of few seconds and it is eKpected that 
the Vanadium Titanium Carbonitrxde particle that formed first must 
have coarsened to some extent in 8 minutes This cannot be 
avoided, however, in this alloy system in which the second phase 
forms by interface controlled mechanism Ct23 The coordinates of 
Figure 3 a(a) are g(p) and p Where p<a/a) is the ratio of the 
edge length to the average edge length and g(p> is the normalized 
distribution function, i a , the fraction of the particles that 
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lie in an interval multiplied by the reciprocal of the interval 
itself 

The particle sire distributions after aging far th, 8h, 27h , 
64h, 216h at 770 C are displayed in Figure 3 8(b), 3 0(c), 3 B(d), 
and 3 0(e) respectively Values of p at which g(jO') is maximum lie 


in the 

range 

of 

1 2 to 2 0 

and the cut off radii ratxo^ 

p , lie 

max 

in the 

range 

of 

1 8ta 2 0 

PSDs after coarsening for 

a 

minuter 

GO 

zr 

64h^ 

137h, 216h and 

343h at 752*^0 ere shown 

in 

Figure 

3 9(a) 

- 3 9(g) 

respectively PSDs for 8 minutes, t h. 

64 

h, 216 


h and 343 h at 733°C are shown in Figure 3 10(a)-3 10(e) 

respectively Values of p at which g(^) is maximum lie in the 

range of 1 0 to 1 8 and their cut off radii ratio p , lie in the 

'^max 

range of I 8-2 O Based on the particle size distribution 

displayed in Figure 3 0 to 3 10 at three different temperature, 
the following observations can be made 

(a) The particle size distributions for most aging condition were 

symmetrical about the mean value, indicating a normal 
distribution rather than a log normal distribution predicted 
by the Lifshiftz, Slyozov and Wagner C2,33 theory The 

coarsening of VC in a Vanadium bearing HSLA Steal is reported 
(4) to have similar PSDs 

(b) The cut off diameter ratio, a/ a is reported to have value 

upto 1 9 instead of the cut off radius ratio p oi 15 as 

predicted by LSW C2,3] theory for the volume diffusion 
controlled coarsening A cut off radius ratio exceeding 1 5 
has also been observed in a number of investigation, namely 
Ni-Co-Al C71, Al-Cu C133 and Fe-C-V C4] 

(c) The peak value of g(p) is observed to fall in the range 



Avg Pttm& Si2e = 7a22 W A\q. t^artide Size = 



Fig 3 8 Particle size distribution, specimen aged at 770°C 










Avg. Particle Size = 242fl6 W A\ft Particle Size = 207^ W 



Fig 3 9 Particle size distribution specimen aged at 

(e) Aged for 137 hours <f) Aged for 216 hours 






Fig 3 Particle -size distribution specimen aged at 733^C 

(b ) Aged far 1 Hour Cc) Aged for 6A- hours 
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1 0-1 4 which IS lower than the theoretically predicted value 
of 2 12 for LSW theory Exception are observed in few casi^s 
where the peak value lies between 1 6 and 1 S 
(d) Instead of being narrow, a number of PSDs were broad near the 
value as displayed in Figure 3 9(c) and 3 10<d) after 
coarsening for a h at 752*^0 and 216 h at 733°C respecti •/'“ly 
The above observation for the coarsening of VTiCN are 
consistent with those reported for vanadium nitride precipitates 
C143 and vanadium carbide C43 Recently Mackay et al 1153 have 
shown that the PSCs during volume diffusion controlled coarsening 
of the }■ phase in Nickel base alloys is a normal rather than a log 
normal distribution 

In an analysis carried out to study the effect of the volume 
fraction on the particle size distribution, Ardell C53 has shown 
that the peak value of decreases with increasing volume 

fraction of the precipitate and the PSD broaden The broadening 
IS symmetrical about value CFig 3 11) At very large volume 

fraction of the precipitate, the PSD for volume diffusion 
controlled coarsening reaches that for interface controlled 
coarsening The PSD, however remains asymmetrical as in the LSW 
C2,33 theory The Lifshiz-Slyozov encounter modified CLSEM3 

theory of Davies et al C73 predicts a normal distribution 
symmetrical about the mean value and is particularly applicable 
for alloy system with a small volume fraction of the precipitate 
CFig 3 121 The basic kinetics equation for the growth remains 

unchanged 

the LSEM theory has any validity then one must observe 


encounter between particles, 


I e , 


coalescence of two or more 




Fiq 3 11 Theoretical PSD modified to accommodate the volume 
fraction of the precipitate CRef 51 



Fig 3 12 Theoretical PSD modified to incorporate 
encounters between particles CRef 73 





neighbouring particles which 


n 


coarsened by receiving salute from 
the matrix During the course ot this investigation we obsenved 
coalescence between particles in a few specimen The encounter 

were between two or more particles in the TEM micrograph of Figure 
3 7(b) In some cases the particles were simply lying one above 
the other Davies et al C73 have shown L shaped and a number of 
rectangular shaped particles of Ni^Al after 44 hours of aging at 
750 C in a Ni~Ca— Al alloy if the LSEM theory of Davies et al 
C73 IS applicable It will, 

(i) Modify the PSD so as to be symmetrical about the mean value 
< 11 ) Increase the cut-off radius ratio to value close to 19 or 
more 


(iii) Broaden the PSD near the mean value, and 

(iv) Lower the peak value of 9(p) 

Although the particle encounter theory of Davies et al C7] 
have some merit consistent with the observation both micro 


structural as well as of particle size distribution, it 
contributes significantly only when the volume fraction of 
vanadium Titanium Carbonitride in HSLA Steels containing 
comparable amount of V, Ti, C and N is reported to be less then 
0 01C163 VTiCN does not nucleate homogeneously but nucleates on 
coherent v/a interface during ferrite transformation by the ledge 
mechanism C17,191 Senerally, a large step (ledge height) is 
followed by several much smaller steps which leads to several 
closely spaced bands of precipitate fallowed by a large space to 
the next band of precipitate C1B3 Thus, a large volume fraction 
of the precipitate is confined within a sheet In pure Fe-V-C 
alloys, an intersect spacing at the order of 200 to 1000 8 is 



reported Cl*?] Because of the large volume fraction of the 
precipitate, encounters between particles within a sheet and 
between particles within a row are expected Prom the shape cf 

the PSDs and micro structural evidence, the encounter modified 
theory of Davies et al C7] deserves some merit 
3 3 Kinetics of Coarsening 

The kinetics equation describing time dependence of the 
particle size is derived for spherical or plate shaped particles 
The Vanadium Titanium Carbonitnde (VTiCN) precipitate is cube 
shaped assumed to have interfacial free energies of the coherent 
face i> The equation derived by Ramakrishna et al C4l is for 
plate shaped particle In their aquation (1 16) ft = i 

( length/ th ickness > is replaced for cube shape particle, the basic 
equations for volume and interface diffusion controlled coarsening 
are identical in form to those of the LSW theory For volume 
diffusion controlled coarsening these are 


H’ ■ H 


„ 0 vv'’ “x" 

8 V e m B 
^ «RT (X^ X") 


(1 20 ) 


Taking the approximate form of the above equation, the 
average value of the particle size, a (S) is plotted against cube 
root of the time (t)^'*^ at the three temperatures used in this 
study Figure 3 13 clearly shows that the average particle size 
data fits the t^^^ proportionality The experimental data for 
volume diffusion controlled growth are shown in Table 3 I 

The coarsening kinetics therefore appears to be volume 
diffusion controlled as also reported for Fe-Mo-C C203, Fe-Ti-Si 
C21] and Fe-V-C C4] alloys The data of Figure 3 13 were fitted 


Tatble 3 1 


Average Particle Size with respect to Time and Temperature 


Temperature 

Ag ing Time 

<h) 

A>'araga Particle Size 

770 

1 

13S 62 

770 

a 

212 02 

770 

27 

254 25 

770 

64 

287 49 

770 

216 

392 40 

752 

t 

108 94 

752 

8 

140 22 

752 

64 

213 93 

752 

137 

1 

242 98 

752 

216 

287 80 

752 

343 

331 31 

733 

1 

70 84 

733 

64 

141 16 

733 

216 

175 90 

733 

343 

231 47 


E 

0 



H6 


to linear least square constant (a^) and slope thus 

determined are shown in Table 3 2 

In an attempt to get some idea of the atomic specie 
controlling the rate of coarsening^ the equation (1 18) was used 
and calculation of the diffusion coefficient D made at the three 

V 

temperatures The constant Y in this equation is dependent upon 
the interfacial energy of the edge with the cn matrix^ the 
parameter a and mole fraction of solute in the at phase in 
equilibrium with a very large size precipitate^ the molar 

volume of the vanadium-titanium carbonitride (VTiCN) and the 
difference (X^ - Xg) 

From the crystal structure (FCC) and the published lattice 

parameters which falls in the range 4 138 to 4 182S for VC^ to 

and in the range 4 3127-4 328S for TiC^ — TiC^ ^ 

Taking the appronimate value of the lattice parameter to be 4 258 

for TiVCN It gives the value of to be 11 558 x 10 ^m^/mole 

For we used B** carbon and determined the mole fraction of 

t? 

carbon in ferrite at various temperatures from the equilibrium 
Fe-C phase diagram This is approximate only The exact value 
can be obtained from the ot/(a+VTiC) boundary in a steel containing 
0 laSVTi and 1 50 Mn Since is very small compared to Xq the 
difference (X^-X^) is taken to be 0 45 for VTiCN The parameter 
a IS treated as -a constant with a value equal to unity ck is not 
^a constant but is dependent upon the particle size as shown by 
Shiflet et al C221 

The interfacial free energy as the coherent broad face is 
expected to be low and comparable to the coherent twin boundary 
energy of 30 mJ/m^ C231 On this basis alone the interfacial 


Table 3 2 


Value of the Constant Ca ) and Slope 

o 


Temperature 

®c 

a 

o 

(8) 

^1/3 

K 

, 3 -1 . 

(ms ) 

770 

78 22 

46 26 

27 5 X 10-^^ 

752 

74 41 

36 67 

13 78 X lO”^^ 




-30 

733 

54 69 

23 4B 

3 6 X 10 









Table 3 3 


Values of K, ard Caleulat-ed 


Temperature 

i 

770®C 

1 . - - - - . 

1 

1 

752°C 

733°C 

if 3 “ 1 

K , ms 

1 27 5xl0“^^ 

13 78x10“^^ 

t 

1 

T i dn-30 

3 6x 1 0 

CO a 

1 486x1 o“^ 

1 

1 72 X lO"^ 

2 135x10“^ 

D , Experimental 

2 “1 
m s 



-18 

1 752x10 

0 830x10“^® 

0 171x10“^® 

r, 2-1 

D , m s 

Ramakr ishana 
at al C43 

— 1 R 

0 16x10 

^0 07x10"'^® 

*•»„ ^_-18 
0 05x10 

I 

r, 2 -1 

D^, m s 1 

Bawen et al C253 

1 

— IS 

282x10 

*168x10“^® 

i .,-.“18 

72x10 

— 

^ 2 1 

D , m s 

Bawen et al C 26 l 

2 7x10"^® 

*1 4x10“'® 

1 

1 

1 

1 

**0 

r, 2-1 

D , tn s 

Mall et al C273 

1 

111x10“^® 

L — j 

67x10"^® 1 

1 

— IS 

12x10 
























2 

energy of the edge can be computed to be about 110 mJ /m The 

2 

average energy of the edge can be taken aa approKimate ly 250 mJ/m 

2 

whxch IS lower than the gram boundary energy (700mJ/m of ferrite 
and interfacial energy of 700mJ/m*" for the Fe/Fe^C interface C24] 
at 710°C 

The value of 0 was calculated from approximate values of the 

V 

parameters in equation (1 IS) as shown in Table 3 J The 

diffusivity values of this investigation are comparable to the 

values obtained by Ramakrishna et al (V in ex— Iron) C43, which is 

very much close to value obtained by Bowen et al <V in ot-Iron) 

C253 and the experimental values differ from those obtained by 

Bowen et al <V in a-Iron) C261 and Moll et al <Ti in cx-Iron) 

C273 by an order of magnitude two The data compare very 

favourably with the experimentally observed values From these 

D values at different temperature, activation energy <Q) and 
v 

can be calculated However, due to insufficient points the 

calculations of and 0^ were not made 



CHAPTER 4 


CONCLUSIONS 


The rata of growth of Titanium Vanadium Carbonitnde in 
HSLA Steel exhibited cube rate law behaviour and x'olurae 
diffusion controlled mechanism predicted by the LSW theory 
of particle coarsening The LSW theory as well as several 
modified theories were ineffective in describing the 
experimentally determined particle size distribution The 
experimentally observed PSDs follow a normal distribution 
rather than a log-normal distribution There is some 
evidence of the modification of the PSDs by encounters 
between particles However, the exact contribution that 
they make in modifying the PSDs is not known The 
experimental particle size distribution differ with respect 
to the peak value, cut off radius ratio and shape of 
predicted PSDs from theoretical considerations Using 
estimated value of the parameter of the kinetics equation, 
the experimental diffusivity, D^, values are comparable with 
the calculated D^ values as they lie within an order of 
magnitude of each other 
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rrax^clers ( rr a __der part i cl a^dens 1 1 y Ck 3 ) nax^d ns part icla^^dfinsity Ck 3 

lower^bord upper_bo d 
k 

uppe _bord + 0 +a g^sizc 

part icl 6 _der 1 1 y EH 0 

) 

avg^ize 00 

tcrCj u )< IP 1 ++) 

C 

a vg 3 z c + dufii 1 C 3 3 /it 
duir rr 3 7 e E 1 3 dum 1 Ci 3 

} 

lowor^tond - 0 0 
upper_Lord 0 ♦ izo 

if (m retr par) 

C 

reiT par rr 
rr u 
> 

) 

numbet^of pack k 1 
tot_p r reir par 


crc t dUftO 0 0 a 0 (<=: 0 rrax_dens)/tot_par 0 8 10 part i cl e.densit y tot.par) 

outr It dtt UMuir pr_of_pacl's tot_par part icle_density ) 

4 , m HtbTOiRAM ***** * ♦**»*♦♦* */' 

prirtrl C > U WANT TO HOU THC 3ZC DISTRIBUTION ON H l^TOCR AH( y/n ) ) 

<iet cf r ( ) 

3 1 { cet ch r 1 ) V i 

C 

rtercjrr) 

piirtri liv THE N rOLO< ( t to s) ? ) 

f, ^f( n ^ilVutT TO UPERTUP) E the ideal LURVE (y^s/noJ 7 ) 

f r / d CT ) 
rr 1- 1 aea r c ( drc 3 ) 

^ \ ^ f le(a 5 SI color 10) 



/ 

H < -I 

/ 


«yt err 
ext t rr 


pTirtft UAIT r ) 

ey t i / ‘^/^^ g/hi to/hg aph </bc: 1 /mV ag/ j np histo >/621 /mkag/f i le out ) 

I ir t f ( ^^L )TTER FILE ( y/n) ? ) 
q t ct 3 ( ) 

1 r ( qe c ( ) y ) 

{ 

q t ch M 

if( IVE THE PLOTTER FILE NAME ) 
ccrt{ I p3 c t_f 1 1 e } 

ftp f ope> (dp w ) 

fr irtf(ftr r plot^file) 

T r 1 t f ( P f 1 7 n ) 
fcl e(f ) 

y ttrr( cat dup >>3 np_pl 1 1 er ) 

VStefT( *3 l/irl' g/histo/hg aph < /s£1 /ini' ag/i np_pl ot ter >/s£1/jnkag/file out ) 
s;sTerr( rn f le ut dup ) 

} 

3 


} 


END OF THE MAIN PROl RAM 

4 I M i * Ji i -{ 


rOUTINi: TO MAKE DATAFILE £0*^ HI TOCRAM 




♦ / 


c t d 1 (xuir xtrax ymin ymax anterval no^points particle^no tot^no^par) 

t)< t xir r xn ymin ym x intfrva) 

ant tc pcints part i ci e_no CdOI tot_no_par 
C 

IL 4 df p tf p n( ) 
j nt 3 

dfr foper ( IZE^DI TRIBUT w ) 
tpr i nt f ( dtp 1 Nr ) 

tpr irt t ( dfp Oi ) 

rp i»Tf(drp f\i/f\t TNtyfNn xirin xuax yrom ymax) 
rprir t(dfp 1 r ) 

fp irtf(dfp a\r no_poi ts) 
for <1 1 i no_point5 i++) 

fen tfCdfr yf r ( O + particle^noCil )/tot_no„par) 
f c I os < df p ) 


H ♦ 


ROUTINE FOR NUMERICAL OUTPUT (STD OUT/TO FILE) »♦♦♦♦♦♦♦*/ 


cutput^delai J (poirt tot_no_par part a cl e_densit y ) 
rt p irt tot_ro_par p rt i cl e_dens i t y CHOI 


t fcairpl e^naire Cl ei 

tic t real_ i eCSOlC 003 overall^avg^size av/g^size 
n t a :rf„ro_par C 03 
fl at lower bond upper^bond 
j r t 3 1 I' n 

t her deci E i < r outf i le C 03 
F 1 L E t > * r )’} v\i ) 

f ( MiITM T IS REDIRECTED TO OME FIUEfy/n) 
ca f( / frde 1 ion) 

J f ( e 1 3 f r ' y ) 
i 

i tf( s n fcarrple_ ame) 
p ircT( A TjCLC. IZ ( Anq Trom)\n\n ) 
r 1 t f ( 
i 

V 1 
for (a 


0 


5 <. no iragrifi ations 1 i++) 



fo 

C 

(] 

1 3 < 1 

iTber_of_p taclesi 3 3 ++) 

P 

f 

c 

t f 

C f t \t 
k f 

r €• l_ 5 i zt Ci3 Cj-1 3 ) 

pri 

f 

} 

f ( 

r ) 


3 

pi 1 rt t ( 

\ ) 


pr 

nt f ( 



. — „ , . 


(* r 1 tt( ^ VERAGF PARTICLE ‘^IZE /8 Sf (Ang trom) \n ov«rftll 
f( unUATF AVEiP t I2£ ( ATISFYINC THE CONDITION ALL PARTICLES SIZE <» 8 0* 
AVERAGE SIZE)«%8 £f (Angstrom)\n avg size \ 

1 i ( d c 11 I n ) ““ ^ 


C 

wifi ive THE OUTPUT FILE NAME ) 
c ft 3 r f 1 I ) 

out f p€*r i L t f i 1 e w ) 
fpr ntf(our s\n air pi e^name ) 
rprirtt((.ut PART 3 CL I 2 C ( An Etron)\n\n ) 

fpr tficut 

1 

for (i 0 i < no ir gnafications 1 i++) 

foi f 1 1 3 < T uiTbor_of_part clasCiI 3++) 

C 

fprintficuf /L Cal t3“ll ) 

If (m n 

C 

tprl tt<oul n ) 

]> 

} 

+ + 


\n\n ) 


} 

f p i rt f i 1 . t 
fpr a rt f ( cuf 
fpr T f ( ou f 
fpr 1 1 ( out 
f p 1 nt 1 ( DU f 
f r r I f ( r 
fpi a I t f ( 0 f 
fra V 
f 3 t f I r u f 
f pi a nt f ( u f 
} 

) 


n ) 


AVERAGE PARTICLE IZE /8 2f (Angstroir) \n overal l_avg_«i 2 » ) 
ULTIMATE AVERAGE <^1ZE /8 df (AngBtrom)\n avg_»i 2 «> 


£ZE KAN E\t\tNUM8ER OF PARTICLE \i ) 


- point - 1 J++) 

t 8 f 8 df\t/a\n O^avg^si za/poa nt ) (i+t)*(2 0*avg_Ba ze/point ) part icle^density Ci+ll ) 


Ut ill ( g oi 0 rclor no^bars) 
av9_ a 

>loi r_b f 


.\n\ 


.\n\ 

An\ 


\n\ 


r 5 i t M j3 

FILE 1 fa n( ) of f a 

t 

pi 1 1 t f f a c 8 £ f avg eazo ) 
afp Irforfinpirito w) 
ipr ti fr IZr^DI TRIBUTXn ) 
tp a f r f p 3 1 0 7\n0 4\ti ) 
fpr a nt f 1 f p \ i j 
f I 0 a < r c_bai 1 a ++ ) 


tp t u r P 

/a r 

colot ) 

Tp r 1 nt r t t p 

^\i ) 


rp r t r a tp 

a Vn 3 \r 

roloi ro)o ) 

tp a r t f r a t p 

i^T ) 


rpr rrfi tp 

/ r j 


fp 1 1 1 1 rp 

1 J 



tpritflifp / f\ ) 

f pi 1 nt T ( 1 f p / r\r ) 

fpriDTfUfp n ) 
f pr i rt f ( i t-p 0 N ) 
fp 1 Tf ( itp 0 \r ) 
tp T f ( 1 f p ( \r ) 
tprirtt(ifp \ i 

f-pr intf ( itp u 0 
fpnrtftirp u y\ ) 
fp3ir(iTp Oul ) 
t pr 3 1 1 1 3 f p \r ^ 

r p 1 I T ( f p 1 o\r ) 
f pr 3 ntf ( If p \n J 

fprintf(ifp o( )\n ) 

tpririt(3fp A<TICL IZZ DJ T'^IBUTION\n ) 

tprirtTtitp r\r ) 

f pr 3 rt f ( 1 f p \r ) 

fpr3rTTiifp 0 3\n ) 

rpri tf^ifp A r ) 

rprirTf( fp 1 x ; 

rp tfufp C\.n } 

fpirtf(fp 07 un) 

rpr3nTfnfp A q a tacle tie /s (A)\n ei 2 «) 

fprlrtfiifp 7\ } 

f pr i nt t ( 1 f p ) l\r ) 
fprirtfi itp hp9 0 \r ) 
f p r 1 nr f ( 3 f p 1 ?\r ) 

fc lofieC ifpj 

ffp foper( 3n|*_plotte w ) 
fprirtffffp IZE^DlSTRI&UTSn ) 

fprtntftfTp 3\r 0 7NnO 4\n ) 
f pr i nt t ( f f p \r ) 

f OT ( a 0 < ro^bars 1 i ++ ) 

f > i { f ( f ft i ID 

f rt tf( ft A ) 

f ri rftrti / i/i\n color colo ) 
f 1\i ) 

frsncfCfj yr) 

fi irtftrr \t ) 

f II If ( f r // f\n ) 

fr 3 nff(rr £ f \n ) 

f nn f(fr> Dn ) 

frstfifr 0 n) 

f r n r r o csn ) 

fj rif f ( f r / I ) 

f ttriflfr Ir ) 

f( rii f t I 0 {) n ) 

f n t f ( f r 0 0 ) 

Prtn(f<fr 0 rt ) 

f rtnrflffi \n 1 

F). ir f (ffi ion) 

f n { n f f) H\n ) 
f} nt tf(fr Q(P)\i ) 

1 tr T(f I PAhTICie IZE DlbTRlBUTtONVi ) 

T j r 1 n f I f r 7 n ) 

r r j t ( f r •} Nt ) 

t T I 0 3\n ) 

f r r t f ( r 4 r ) 
f r 3 1 M ( f r 1 n ) 
t 3 r t f ) 

f 1 t f ( f r 17 Nn ) 

ijrfir Aq Patcle tr 

f 3 Iff 7Kr > 

f r jr f I r Hi ) 

cl f 


yis (A)\n lie) 


m ^ e_ic.e leu c ( deci ) 
thfir dec I 3 
( 

FILE 1 f j -if nC ) 
r I t V 

iff roFe^< / 1 nl'a9/Ka t o/ 1 DE AL_CUR VE w ) 

f r I ( r \ d 1 ) 

f r i { r ( f t 7\n ) 

f f 0 0 ^1*^71 x + 001) 

t f ( I ) & f ( X < 0 039^^9) II (X > 0 14011))) 


( 0 t 87 0 61 0103 + x 3 0S45 1 4G&49’**xpt3 0 l 09 ix))-S 36193 •xp(4 0*lo&tx)) 

y 7 ^ O’^j expi 0-flOQ<x)> 0 S 3 4 «xp(6 0 lc»g<x))*f0 80£699^«Xp(7 Q loglxll^S 27141 *Xp(8 O^loyCxli 

fp r i f t 3 rp f T/ f Sn x y ) 
j 

el if ( (y > 1 ) ^4 (X < 1 37)) 

C 

y 14 f ( 3 6727*x-^ 6)9 fX x + l£ 6906*expf3 0 ]oc1x))-13 7876*#xp<4 0 1 oc ( x ) ) 

I 4“* -s? exp< 0 loq(x)) 1 018^ €fXp(6 0^ 1 09 ( x ) ) + 0 087391 txpIT 0 >»lo 9 (x)) 9 866£9**xpt8 0 lo^lx)) 

4o7 exp(9 0>1cvq(x)) 
fpr 1 Tffirp /f t/r\n x y) 


tprirtffafp 1 
f p r 3 nt r ( f p 1 

TprirtTiiTp 1 

f print f i i tp 1 
f c 1 Gse i 1 f p 3 


1 14\t 0 017 \n ) 
4 \t 0 01 j\n ) 
4l7\t0 00 \n ) 
434 \t 0 0000\n ) 



